The Bethe ansatz as a matrix product ansatz 
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The Bethe ansatz in its several formulations is the common tool for the exact solution of one di- 
mensional quantum Hamiltonians. This ansatz asserts that the several eigenfunctions of the Hamil- 
tonians are given in terms of a sum of permutations of plane waves. We present results that induce us 
to expect that, alternatively, the eigenfunctions of all the exact integrable quantum chains can also 
be expressed by a matrix product ansatz. In this ansatz the several components of the eigenfunctions 
are obtained through the algebraic properties of properly defined matrices. This ansatz allows an 
unified formulation of several exact integrable Hamiltonians. We show how to formulate this ansatz 
for a huge family of quantum chains like the anisotropic Heisenberg model, Fateev-Zamolodchikov 
model, Izergin-Korepin model, Sutherland model, t — J model, Hubbard model, etc. 



Since the pioneering work of Bethe in 1931 [1] the 
Bethe ansatz and its generalizations proved to be a quite 
efficient tool in the description of the eigenvectors of a 
huge variety of one-dimensional quantum chains and two 
dimensional transfer matrices (see e. g. [2]- [5] for re- 
views). On the other hand in the last two decades [6]- 
[9] it has been shown that some special quantum chains, 
although not integrable through the Bethe ansatz have 
the components of its ground-state wavefunctions given 
in terms of a product of matrices. In this matrix prod- 
uct ansatz (MPA), apart from a normalization constant, 
these components are fixed by the algebraic properties 
of the matrices defining the MPA. In another context a 
MPA has also been applied quite successfully to the eval- 
uation of the stationary distribution of probabilities of 
some stochastic models in one dimension [10]. The time 
fluctuations of these stochastic models are described by 
the ground state wavefunction of a related spin Hamilto- 
nian. The simplest example is the one-dimensional asym- 
metric exclusion process [10] whose related spin Hamilto- 
nian is the anisotropic Heisenberg chain, or XXZ chain, 
with appropriate boundary fields [11]. The stationary 
properties of the model are given in terms of the alge- 
braic relations of the matrices appearing on the MPA. 
This ansatz was used in a variety of problems including 
interface growth [12], boundary induced phase transitions 
[10], [13]- [15], the dynamics of shocks [16] or traffic flow 
[17]. 

An important development of the MPA that appeared 
in the context of stochastic models is the dynamical ma- 
trix product ansatz (DMPK) [18,30]. This ansatz allows, 
whenever it is valid, the calculation of the probability 
densities, of the stochastic system, at arbitrary times. In 
the related spin Hamiltonian this DMPA asserts that not 
only the ground-state wave function, as in the standard 
MPA, but an arbitrary wavefunction have components 
expressed in terms of a matrix product ansatz whose ma- 
trices, in distinction of the standard MPA, are now time 
dependent. This DMPA was shown originally [18], [19] to 



be valid for the problem of asymmetric diffusion of par- 
ticles on the lattice. More recently [20] it was also shown 
that this DMPA can also be formulated in the context of 
stochastic models with two species of hard-core particles. 
The condition of validity of the DMPA reproduces the 
subspace of parameters where the model is known to be 
exactly integrable. This fact induce us to conjecture that 
all the Hamiltonians, related or not to stochastic mod- 
els, that are solvable through the Bethe ansatz may also 
be solvable by an appropriate MPA. This would mean 
that the components of the eigenfunctions of the exact 
integrable models, that according to the Bethe ansatz, 
are normally given by a combination of plane waves, can 
also be obtained from the algebraic properties of the ma- 
trices defining the MPA. 

In this Letter we are going to show how the eigen- 
spectra of the exact integrable quantum chains, with 
global conservation laws, can be obtained by an appro- 
priate matrix-product ansatz. In this way we were able 
to obtain the integrability of several well known exact 
integrable models. Among the models with one global 
conservation law we have the XXZ chain [21], the spin- 
S Fateev-Zamolodchikov model [22], the Izergin-Korepin 
model [23], the solvable spin-1 model of Ref. [24], etc, 
and among the models with two global conservations we 
have the supersymmetric t — J model [25], the spin-1 
Sutherland and Perk-Schultz models [26,27], the Hub- 
bard model [28], as well the two-parameter integrable 
model presented in [29]. Moreover the matrix-product 
ansatz we propose enable us to show, with little effort, 
how to extend the above mentioned models, by including 
arbitrary range of hard-core interactions, without losing 
their exact integrability. 

For brevity and in order to illustrate the proposed 
MPA we are going to present here two examples. The 
solution of the XXZ chain with arbitrary hard-core in- 
teractions among the up spins, as an example of a model 
with one global conservation law, and the solution of the 
Hubbard model as an example of a model with two global 



1 



conservations laws. 

The Haniiltonian of the XXZ chain with a hard-core 
exclusion of S sites (5* = 1,2,...) is the anisotropic 
Heisenberg chain (anisotropy A) where any two up spins 
are not allowed to occupy lattice sites at distances smaller 
than S. This Hamiltonian in a L-sites periodic chain is 
given by 



Hs = -Vs + ^i^i+i + ^'^M+i)'Ps, (1) 

i=l 

where a^^a^^a^ are spin-| Pauli matrices and Vs is a 
projector that projects out the configurations where any 
two up spins are at distances smaller than S. The case 
S = \ corresponds to the standard exactly solvable XXZ 
chain [21]. The conserved charge associated to the global 
conservation law of (1) is the z- component of the total 
magnetization, or equivalently the number of up spins 
(particles). The translation invariance of the lattice also 
ensures that the momentum is also a good quantum num- 
ber. 

The ansatz we propose states that any of the wave- 
functions IV'n.p > in the sector with n spins up (n = 
0, 1,2, . . .) and momentum P ^ = 0, 1, 2, . . . , i- 1) 
is given by a matrix product ansatz, i. e., their ampli- 
tudes are given by the trace of the matrix product: 



\tPn,p>= Tr(£^i-iA£;^^-^i-iA-. 

Xi,...,Xrt 



U£^-^"np)|xi,...,x„ >, 



(2) 



where . . . ,a;„ > denotes the configurations with up 
spins at (xi, . . . , Xn) and the symbol (*) in the sum means 
the restriction to the configurations where \xi^i—Xi\ > S. 
The algebraic properties of the matrices A, E and ftp will 
be fixed by the eigenvalue equation 



(3) 



The algebraic relations of the matrices A and E with 
rip will fix the momentum of the eigenfunction. The 
fact that IV'n.p > has a momentum P imply that the 
ratio of the amplitudes corresponding to configurations 

\xi, . . . ,Xn > and + 1, . . . , a;„ + 1 > is e~*^ and con- 
sequently from (2) we have the following relations 



AQp = e-'^QpA, EQp = e-'^QpA. 



(4) 



The eigenvahic equation (3) when applied to the compo- 
nents of \ipn,p > where all the up spins are at distances 
larger than S give us the constraint 

e„Tr(- • • . . . AE^-^'^np) 



A convenient solution of this last equation is obtained by 
identifying the matrices A as composed by n spectral- 
parameter dependent matrices 



(6) 



where the matrices A/^^ obey the commutation relations 

EAk, = e^'^^Ak^E, Ak ^p = e-^P^'^-^HlpAk^, (7) 

and kj {j = 1, . . . ,n) are in general complex numbers 
unknown a priori [31]. The energy and momentum are 
given respectively by 

n n 
j=l i=l 

The eigenvalue equation when applied to the other 

components of \ijjn,p > will give the commutation re- 
lations of the matrices {^/c^} among themselves. In 
fact this algebra is obtained from the components where 

any pair of up spins are located at the closest positions 
("matching conditions") Xj and xj-^-i = Xj + S, namely 

Ak^Ak, = 5(fc„ ki)Ak^^Ak, {i + j), A\ = 0, 

giikj+h) _|_ X _ 2Ae*'°^ 
S{kj, h) = - ^.^^^^^^^ ^ ^ _ ^^^.^^ ,i,j = l,...,n. (9) 

No new algebraic relations appear for the matrices Ak^ 
and the associativity of the algebra (7) and (9) follows 
from property S{ki, kj)S{kj, ki) = 1. The cyclic property 
of the trace in (2), together with the algebraic relations 
(9) give us, for each kj (j =1,. . . ,n), the conditions 

Tr(Afc, ■ • • Afe„^^-"('^-i)Op) = e-''=i(^-"'5+")e-'^(^-i) 

n 

x([] S{k,, fcO)Tr(Afe, ■ • • Afe„£;^-"(^-i)np), 
1=1 

that fix the up to now free complex spectral parameters 



ikjL 



n5(fc,-,fciK^'^-''^^^-'\ (i = l,...,n). (10) 



(=1 



The wavefunctions, obtained by using (2) (8) and (9) can 
be written as a combination of plane waves ("wave num- 
bers" {kj}), and apart from a normalization constant co- 
incides with those obtained by th(> standard Bethe ansatz 
[21] for 5 = 1 or arbitrary values of S [32]. 

As a second example we consider the standard Hub- 
bard model whose Hamiltonian in a periodic lattice with 
L sites is given by 



= - E[Tr(- • • • • • AE^-''-np) Hh = ^(c+^Cj+i,<, + c+^^^„Cj,„) + uY,rij^-r]j,+ , 



i=l 



+Tr(- • • £:^i-^*-M^^i+i-^i-2^ . . . AE^'^'^np)]. (5) 



(11) 
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where c^^ are creation operators of electrons of spin 

cr = ± at site j, and rij^a = "-/o-^jio- are the number 
operators of electrons of spin a at the site j. In this case 
we have two global conservation laws n±, corresponding 
to the number of electrons with spin cr = ±, and the toal 
number of electrons is n = n+ + n_ . In order to form 
our MPA we associate, as before, the matrix E to the 
empty sites, the matrices X~ to the single occupied 
sites with electrons with spin up and down, respectively, 
and the matrix X° = X~^E~^X~ to the sites with dou- 
ble occupancy. Our MPA asserts that the components 
of the cigcnfunction \4'n,P > of energy e„ and momen- 
tum P = 2itI/L [l = 0, 1, . . . , L — 1) corresponding to 
the configuration \x\,Q\] . . . ,Xn,Qn > where the non- 
empty sites (.xi, . . . , .T„) have occupation (Qi, . . . , Qn) 
{Qi = +1 ^^ 0) will be given by the trace 

Tr(i;^i-^X'5i£:^2-a;i-i^Q2 . . . x^^E^-'-^^^lp). (12) 

The momentum of the state P, as in (4), is fixed by im- 
posing the commutation relations X'^Vlp = e^^^ilpX^ , 
[Q = +,—,0) and Eflp = e~^^^pE. The eigenvalue 
equation Hh\'^n,p >= ^n\'^n,p > will provide the alge- 
braic relations of the matrices E and X'^. 

The components of the wavcfunction corresponding to 
the configurations where all the particles are at distances 
\xi+i —Xi\ > 1 will give a generalization of (2) whose solu- 
tion is obtained by introducing the convenient n spectral- 
parameter dependent matrices 

X« = f^EXg, {Q = +,-), (13) 

whose commutation relations with the matrices E and 

flp are 

EX^^ = e^'^^X^.E, X^^np = npX^^. (14) 

The energy and momentum in terms of these unknown 
complex spectral parameters are given, as in (8), by 

n n 

The components where the particles occupy the clos- 
est positions and those where we have double occupancy, 
give us, by using (14) and (15) the algebraic relations 

= s22iki, kj)x^xl {X^f = 0, 

where Q' = —Q and Q = ±. The relations (14) and (16) 
define completely the algebra whose structural constants 



arc the well know S'-matrix of the Hubbard model [28]. 
Since the several components of the wavcfunction should 
be uniquely related the above algebra should be associa- 
tive. This associativity imply that the above S'-matrix 
should satisfy the Yang-Baxter relations [33,2], which is 
indeed the case [28]. The components of the wavcfunc- 
tion corresponding to the configurations where we have 
three or four particles in next-neigboring sites would give 
in principle new relations involving three or four opera- 
tors X^. . These new relations arc however consequences 
of the above relations (14) and (16). The cychc property 
of the trace in (12) and the algebraic relations (14) and 
(16) will imply 

E n{Qh{Q'm{x^}---x^}E^np). (17) 

Q[,-,Q'„ 
where 

n 

Q"i,...,Q"„ i=l 

is the transfer matrix of a non-homogeneous six-vertex 
model defined on a cylinder of perimeter n and with 
Boltzman weights given by the S-matrices defined 
in (16). The eigenvalues of this auxiliary problem 
^n{ki;ki, . . . , kn) can be obtained in a standard way by 
the coordinate-Bethe ansatz [21] or by the quantum in- 
verse scattering method [34]. Using these eigenvalues in 
the relation (17) the spectral parameters {kj} will be 
fixed by the solutions of the system of equations 

e-ifc.i = (_l)"A„(fcj;fci,...,fc„), (j = l,...,n). (19) 

These last equations coincide with the Bethe-j4nsate 
equations derived through the standard coordinate Bethe 
ansatz [28]. 

Generalizations of our MPA is quite simple to imple- 
ment [35] . For example the solution of the excluded vol- 
ume Hubbard model where electrons with spin up (down) 
exclude other electrons at S+ (S_) sites on its right, but 
allows double occupancy at any site, is just obtained by 
changing in (13) X^ = J2j E^'^^'^X^^ (Q = ±)- 

We have also shown [35] that the above MPA also 
works for the other known exact integrable models with 
two conserved global quantities (C/(l) (8i U{1)), like the 
Essler-Korepin-Schoutens model [36] or the generalized 
two-parameter integrable model introduced in [29]. In 
those last cases the same MPA presented above for the 
Hubbard model apply except that now the matrices X'^ 
associated to the sites with double occupancy are given 
in terms of new operators X^ = Y~^E~^Y~ , and we were 
able to rederive all the results obtained for those mod- 
els through the standard Bethe ansatz. Hamiltonians 
with two global conservations that do not allow double 
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occupancy like the stochastic Hamiltonian associated to 
the asymmetric diffusion with two Icinds of particles [32] , 
or the supersymmetric i — J quantum chain [25] or the 
SU{3) Sutherland and Perk-Schultz models [26,27], are 
obtained through the MPA (12), without the matrices 

Models of spin one with a single global conservation 
law (conservation of the 2;-magnetization) like the Fateev- 
Zamolodchikov model [22], the Izcrgin-Korcpin model 
[23], or the solvable spin-1 Hamiltonian introduced in 
[29], have their solutions given by a MPA similar as in (2), 
where now we associate the matrices E, A and BE~^B 
to the sites occupied by particles with spin S^- — ~1, 
or +1, respectively. 

We have also obtained an appropriate extension of the 
presented MPA to non periodic, but exact integrable 
boundaries as, for example, the XXZ chain with surface 
fields. 

It is interesting to notice that in the cases of exact inte- 
grable Hamiltonians associated to stochastic models, as 
in [11,30], since we can write all eigenfunctions in a MPA, 
our results imply that we can equivalently write at any 
time the probability distribution of the model in terms of 
a time-dependent MPA, as happens in the DMPA [18]. 

In conclusion, we have shown that the eigenfunctions 
of a huge variety of exact integrable quantum chains can 
be expressed on an unified way in terms of a matrix prod- 
uct ansatz, whose; matriccis satisfy an associative algebra. 
The associativity of the algebra that warranties the exact 
integrability of the model is a consequence of the Yang- 
Baxter relations. 

This work has been partially supported by FAPESP, 
CNPq and CAPES (Brazilian agencies). 
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